CHAPTER 7 


THEORIES OF BEARING CAPACITY 
AND SETTLEMENT 


7.1 INTRODUCTION 

Testing of soil in the laboratory, and perhaps in the held, to obtain properties 
is required to allow the computation of bearing capacity and settlement. The 
reader is referred to Chapter 3 for details on the determination of the required 
parameters of the soil at a site. 

The problems to be addressed in this chapter can be illustrated by the 
shallow foundation of given lateral dimensions resting on soil, as shown in 
Figure 7.1a. The first of two problems facing the engineer is to find the unit 
vertical load on a shallow foundation, or on the base of a deep foundation, 
that will cause the foundation to settle precipitously or to collapse. The unit 
load at failure is termed the bearing capacity. With regard to bearing capacity, 
the following comments are pertinent: 

1. Any deformation of the foundation itself is negligible or disregarded. 

2. The stress-strain curve for the soil is as shown in Figure 7.1b. 

3. The base of the foundation may be smooth or rough. 

4. The soil is homogeneous through the semi-infinite region along and 
below the foundation. 

5. The loading is increased slowly with no vibration. 

6. There is no interaction with nearby foundations. In spite of the con¬ 
straints on computing the values, the concept of bearing capacity has 
been used for the design of foundations for decades and presently re¬ 
mains in extensive use. 
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Figure 7.1 Example of footing and stress-strain curve used in developing bearing- 
capacity equations. 


Investigators have noted that the bearing capacity equations will vary if 
the failure is general where symmetrical failure surfaces develop below the 
base of the foundation; is local where a failure occurs due to excessive set¬ 
tlement of the foundation; or is a punch-through where the foundation 
punches through a strong surface layer and causes the weak soil below to 
fail. The equations for general shear failure are most important. 

The second of the two problems addressed here is the computation of 
settlement of a foundation such as that shown in Figure 7.1. Two types of 
settlement are noted: immediate or short-term settlement and long-term set¬ 
tlement due to the consolidation of saturated clays. The immediate settlement 
of foundations on sands of loose or medium density is relatively so large that 
settlement controls, termed a local shear failure, as noted above, and a general 
bearing capacity failure does not occur. Immediate settlement of foundations 
on sand and clay is discussed in Chapter 9. 

Equations and procedures for dealing with long-term settlement due to 
consolidation are presented here. The settlement of deep foundations is dis¬ 
cussed in chapters that deal specifically with piles and drilled shafts. 

The finite-element method (FEM) discussed in Chapter 5 provides valuable 
information to the engineer on both bearing capacity and settlement, as dem¬ 
onstrated by the example solution presented. FEM can now be implemented 
on most personal computers, rather than on large mainframes, and will play 
a much greater role in geotechnical engineering as methods for modeling the 
behavior of soil are perfected. Leshchinsky and Marcozzi (1990) performed 
small-scale experiments with flexible and rigid footings and noted that the 
flexible footings performed better than the rigid ones. Rather than use the 
performance of expensive tests with full-sized footings, FEM can be used to 
study the flexibility of footings. 
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7.2 TERZAGHI’S EQUATIONS FOR BEARING CAPACITY 

The equations developed by Terzaghi (1943) have been in use for a long time 
and continue to be used by some engineers. His development made use of 
studies by Rankine (1857), Prandtl (1920), and Reissner (1924). The two- 
dimensional model employed by Terzaghi (Figure 7.2) was a long strip foot¬ 
ing with a width of 2 B (later, authors simply use B for the width of footing) 
and a depth of the base of the footing below the ground surface of D t . The 
base of the footing was either smooth or rough, with the shape of the solution 
for a rough footing shown in the figure. The wedge a, b, d is assumed to 
move down with the footing, the surface d, e is defined by a log spiral and 
the soil along that surface is assumed to have a failure in shear, and the zone 
a, e, f is assumed to be in a plastic state defined by the Rankine equation. 

The equations developed by Terzaghi for bearing capacity factors are as 
follows: 

Q d = 2B(cN r + yD f N q + yBN y ) (general shear failure) (7.1) 

Q' d = 2B(% cN' c + yBN' y + yBN’ y ) (local shear failure) (7.2) 

where 

Q d = load at failure of the strip footing, 
c = cohesion intercept, 
y = unit weight of soil, and 

N c , N' c , N q , N' q , N r N' y = bearing capacity factors shown graphically in 

Figure 7.3. 

The differences in settlement for a general shear failure and a local shear 
failure are shown in Figure 7.4, where the general shear failure is depicted 
by the solid line. The footing on loose soil is expected to settle a large amount 
compared to the footing on dense soil, which Terzaghi elected to reflect by 



Figure 7.2 Model used by Terzaghi in developing valves for bearing-capacity factors 
(from Terzaghi, 1943, p. 121). 
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Figure 7.3 Terzaghi’s bearing-capacity factors (from Terzaghi, 1943, p. 125). 


reducing the values of the bearing capacity factors for loose soil. A proposed 
reduction in bearing capacity to deal with excessive settlement is disregarded 
here because Chapter 9 deals with short-term settlement of shallow founda¬ 
tions. The engineer may reduce the load on the foundation if the immediate 
settlement is deemed excessive. 


7.3 REVISED EQUATIONS FOR BEARING CAPACITY 

A number of authors have made proposals for bearing capacity factors, in¬ 
cluding Caquot and Kerisel (1953), Meyerhof (1963), Hansen and Christensen 
(1969), Hansen (1970), and Vesic (1973). The basic form of the bearing- 
capacity equation (Eq. 7.1), proposed by Terzaghi, has been accepted by most 
subsequent investigators; however, two modifications have been suggested: 
(1) improved analysis of the model proposed by Terzaghi (Figure 7.2), and 


Unit Load 



Figure 7.4 Assumed curves for unit load versus settlement for a footing on dense 
soil, Q and on loose soil, C 2 (from Terzaghi, 1943, p. 118). 
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(2) an extension of the method to include a number of factors such as inclined 
loading and a rectangular footing rather than a strip. 

In addition, a number of other studies have been undertaken, such as those 
using FEM, to study the upper- and lower-bound values of the bearing ca¬ 
pacity (Ukritchon et al., 2003). Hansen’s method has been selected and is 
presented in detail because it is comprehensive and because the analytical 
results agree reasonably well with experimental results. 


7.4 EXTENDED FORMULAS FOR BEARING CAPACITY 
BY J. BRINCH HANSEN 

Hansen employed the same basic equation as Terzaghi except that the width 
of the footing is B instead of 2 B as employed by Terzaghi. Hansen’s basic 
equation is 

Y = \ ^ BN y + y°f N , + cN < ( 7 - 3 ) 

where 

Q d = ultimate bearing capacity, 

B = width of foundation, and 

y = unit weight of soil (use y' for submerged unit weight). 

The equations below are for the bearing-capacity factors N q and N c as ex¬ 
pressed by exact formulas developed by Prandtl (1920): 

N q = e"*** tan 2 ^45 + (7.4) 

N c = (N q - l)cot 4> (7.5) 

(Note: if 4> = 0, N c = it + 2.) 

Hansen and Christensen (1969) presented a graph for N y as a function of 
the friction angle, S, between the base of the footing and the sand (Figure 
7.5). If the base is rough, S = 4>, and the values can be read from the indicated 
curve in Figure 7.5. Footings with a rough base, as would occur with concrete 
poured on a base of sand, are usually assumed in design. One of the authors 
witnessed construction in Moscow where concrete footings were in a plant 
and were being trucked to the job site. Extremely low temperatures for many 
months prevent casting of concrete in the held. If a designer wishes to account 
for low friction between the base of a footing and the sand, the curves in 
Figure 7.5 may be used, with the values for 8 = 0 indicating a perfectly 
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Figure 7.5 Bearing-capacity factor N found for strip footings as a function of fric¬ 
tion angles <p for sand and 8 for subsurface of footing (from Hansen and Christensen, 
1969). 


smooth footing. The following equation is for a perfectly rough footing and 
is commonly used in design: 

N y = 1.5 (N q - l)cot $ (7.6) 

Hansen’s bearing-capacity factors for the strip footing are presented in 
graphical form in Figure 7.6 and tabulated in Table 7.1. Perhaps Hansen’s 
most important contribution was to extend the equations for bearing capacity 



Bearing Capacity Factor 

Figure 7.6 Hansen’s bearing-capacity factors. 
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TABLE 7.1 Hansen’s Bearing-Capacity Factors 


<p (degrees) 


N c 

N y 

0 

1.0 

1.5 

0.0 

5 

1.6 

6.5 

0.1 

10 

2.5 

8.3 

0.4 

15 

3.9 

11.0 

1.2 

20 

6.4 

14.8 

2.9 

25 

10.7 

20.7 

6.8 

30 

18.4 

30.1 

15.1 

31 

20.6 

32.7 

17.7 

32 

23.2 

35.5 

20.8 

33 

26.1 

38.6 

24.4 

34 

29.4 

42.2 

28.8 

35 

33.3 

46.1 

33.9 

36 

37.8 

50.6 

40.1 

37 

42.9 

55.6 

47.4 

38 

48.9 

61.4 

56.2 

39 

56.0 

67.9 

66.8 

40 

64.2 

75.3 

79.5 

42 

85.4 

93.7 

114.0 

44 

115.3 

118.4 

165.6 

46 

158.5 

152.1 

244.6 

48 

222.3 

199.3 

368.7 

50 

319.1 

266.9 

568.6 


to deal with deviations from the simple case of a strip footing. The load may 
be eccentric, inclined, or both. The base of the foundation is usually placed 
at depth D f below the ground surface. The foundation always has a limited 
length L and its shape may not be rectangular. Finally, the base of the foun¬ 
dation and the ground surface may be inclined. The equation for Hansen’s 


extended formula is 


-f = n y BN y s ydyiy b ygy + yD f N q s q d q i q b q g q + cN c s c d c i c b c g c (7.7) 


where 


A = area of footing, 

5 = shape factors, 
d = depth factors, 
i = inclination factors, 
b = base inclination factors, and 
g = ground inclination factors. 
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In the special case where $ = 0, Hansen writes that it is theoretically more 
correct to introduce additive factors, and Eq. 7.8 presents the result: 

^-' = (77 + 2K(1 + X + d“ - i“ - K - & (7.8) 

where s u = undrained shear strength of the clay. The definitions of the mod¬ 
ification factors are shown in the following sections. Hansen stated that when 
the modifications occur one at a time, a simple analytical solution or results 
from experiments can be used; however, when all of the factors are used 
together for more complicated cases, the resulting computation will be an 
approximation. 

All of the loads acting above the base of the foundation are to be combined 
into one resultant, with a vertical component, V, acting normal to the base 
and a horizontal component H acting in the base. To account for eccentric 
loading, the foundation is re-configured so that the resultant intersects the 
base at a point called the load center. If the foundation has an irregular shape, 
a rectangular foundation meeting the above conditions is employed. 


7.4.1 Eccentricity 

With regard to eccentric loading on a shallow foundation, Hansen recom¬ 
mended that the foundation be reconfigured so that the eccentricity is elimi¬ 
nated on the foundation to be analyzed. The scheme proposed by Hansen is 
shown in Figure 7.7, where the original foundation is indicated by the solid 
lines. The point of application of the loading on the base of the foundation 



Figure 7.7 Sketch showing reconfiguration of the shape of a foundation to account 
for eccentric loading. 
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is shown by the dot, and the foundation to be analyzed is shown by the dashed 
lines. The eccentricity is exaggerated for purposes of illustration. 

For foundations of other shapes with eccentric loads, the best possible 
rectangular foundation is configured with no eccentricity to replace the par¬ 
ticular shape. 


7.4.2 Load Inclination Factors 

An inclination of the load will always mean a reduced bearing capacity, and 
solutions can be developed by revising the model shown in Figure 7.2. The 
following simple empirical formulas were proposed by Hansen: 

lq= l ~ °' 5 (V + Ac cot 0) 

• = h 07 H 
ly _ (V + Ac cot 4>) 

For the case where 4> = 0, the following equation may be used: 
i a c = 0.5 - 0.5 ^1 

where 


H 

Ay" 


(7.11) 


(7.9) 

(7.10) 


H = component of the load parallel to the base and 

V = component of the load perpendicular to the base. 

Equations 7.9 and 7.10 may not be used if the quantity inside the brackets 
becomes negative. 

The above equations are valid for a horizontal force where H is equal to 
H B that acts parallel to the short sides B of the equivalent effective rectangle. 
If H B is substituted for //. the computed factors may be termed i qB , i yB , and 
i a cB . In the more general case, a force component H, also exists, acting parallel 
with the long sides L. Another set of factors may be computed by substituting 
H l for H, and factors will result that may be termed i qL , i yL , and i a cL . The first 
set of factors with the second subscript B may be used to investigate possible 
failure along the long sides L, and the second set of factors with the second 
subscript L may be used to investigate possible failure along the short sides 
B. One of the sets of analyses will control the amount of load that may be 
sustained by the foundation. 
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7.4.3 Base and Ground Inclination 


The sketches in Figure 7.8 defined the inclination v of the base of the foun¬ 
dation and the inclination (3 of the ground. In the case of 4> — 0, the following 
equations were found: 

, 2v v 

b“ = -=- 

c tt + 2 147 

2/3 j3 

e a = — = -2— 

Sc tt + 2 147 

where v, /3, and 147 are in degrees. 

For values of 4> other than zero, the following equation applies: 

b q = e - 2vtim<l> (7.14) 

Equations 7.11 through 7.13 may be used for only positive values of v and 
/3, with /3 smaller than 4>, and v plus (3 must not exceed 90°. 


(7.12) 

(7.13) 


7.4.4 Shape Factors 

Hansen recommended the following equations to deal with footings with rec¬ 
tangular shapes: 



V 



Figure 7.8 Sketches defining inclinations (a) of base v and (b) of ground /3 (from 
Hansen, 1970). 
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s y = 1 - 0.4 ^ 

(7.15) 

s q = 1 + sin d> - 

(7.16) 

s “ c = 0-2 z 

(7.17) 


The above shape factors are only for vertical loads. For inclined loads, 
formulas must be introduced to account for the inclination. Failure may occur 
along the long sides or along the short sides; thus, two sets of equations are 
required. 


J qB 




3 7 B 


>yL 


= 0.2 i“ cB | 

(7.18) 

L 

= o 9 i a — 

■ CL B 

(7.19) 

= 1 + sin 4> 

(7.20) 

= 1 + sin 4> 

B 

(7.21) 

= 1-0.4 ^ 

Tj-yZ. 

(7.22) 

LLr 

= 1-0.4-^ 

B 

(7.23) 


For Eqs. 7.22 and 7.23, the special rule must be followed that a value 
exceeding 0.6 must always be used. 

7.4.5 Depth Effect 

Foundations are always placed at depth D f below the ground surface, so the 
term yl), is used to indicate the soil above the base of the foundation. The 
shear strength of the overburden soil is assumed to be the same as that of the 
soil below the base of the footing. If the overburden soil is weaker than 
the foundation soil, the depth effect may be reduced or ignored entirely. 


7.4.6 Depth Factors 

The depth factor for the effect of the weight of the soil is always equal to 
unity: 
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The result in Eq. 7.29 agrees with the well-known value found for the 
resistance of the tip for piles founded in clay. The following formula applies 
for other values of friction angles: 

% = yD f N (1 + sin <£)[1 + 7r tan <p(l - sin 4>) 2 ] (7.30) 

A v 

Hansen noted that for friction angles between 30° and 40° the above equa¬ 
tion yields results in good agreement with the Danish experience for point 
resistance for piles in sand, provided that d> is taken as the friction angle in 
plane strain. 

For the usual case of failure along the long sides L on the base, Eqs. 7.29 
and 7.30 may be used with Eqs. 7.25 and 7.26, yielding values of the depth 
factors d“. B and d qB . For the investigation of a possible failure along the short 
sides, 5, equations for another set of depth factors must be used: 
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d a cL = 0.4 arctan -j- (7.31) 

„ D, 

d qL = 1+2 tan <p( 1 - sin 4>) 2 arctan -j- (7.32) 

7.4.7 General Formulas 

In many instances, the horizontal force will have a component H B parallel 
with the short sides B and a component H L parallel with the long sides L, 
and the following formulas must be used: 

\ y N y Bs W 7 B b -y + (y° f + c cot 4>)N q d qR s qH i qH b q - c cot c/> (7.33a) 

or 

\ y N y Bs yd y r b y + (7°/ + c cot 4>)N, / qL s qL i qL t> q - c cot (f> (7.33b) 

The above equations may be used as follows. Of the two possibilities for 
the terms including y, the upper term should be used when Bi yB < Li y , and 
the lower term should be used when Bi y[i > Li L . A check on the proper 
choice is that s y > 0.6. Of the two possibilities for the terms including yD f , 
the one giving the smallest value must always be chosen. 

In the special case where cf> = 0, the smallest of the values from the 
following two equations may be used: 

^-'<(77 + 2)c„(l + s“ cB + d“ B - i‘‘ B - b‘‘ - g a c ) (7.34a) 
or 

^-'<(77 + 2)c„(l + S‘‘ L + d a cL - i a cL - K - g!) (7.34b) 

7.4.8 Passive Earth Pressure 

If sufficient lateral movements of a footing occur to develop passive earth 
pressure on the edge of the footing, the bearing-capacity equations will be 
affected. Hansen’s studies showed that the movements necessary to develop 
bearing capacity and passive earth pressure are of the same order of magni¬ 
tude. The movement associated with bearing capacity is discussed in Chapter 
9. In special cases, the movements required to develop bearing capacity and 
passive earth pressure can be computed by FEM. 
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With respect to design, for mainly vertical loads no passive pressure is 
employed. For mainly horizontal loads, passive pressure must be employed; 
however, the details of construction play an important role. If an oversized 
excavation is made for the foundation, the nearby soil may have been dis¬ 
turbed and/or the backfill may have been placed poorly so that passive earth 
pressure does not develop. 

7.4.9 Soil Parameters 

In the case of saturated clays, the relevant value c u of undrained shear strength 
must be employed for short-term behavior; for long-term behavior, the param¬ 
eters (b and c must be employed as determined from drained triaxial tests. 

For sand, c can be assumed to equal zero. The equations shown in the 
sections above are based on the friction angle in plane strain. The following 
identity is recommended: 

= 11 <t>, r (7-35) 

where 

4> pl = friction angle for plane strain conditions, and 
4> tr = friction angle determined from triaxial tests. 

7.4.10 Example Computations 

Examples from Hansen (1970) The results from two tests of a foundation 
on sand were reported by Muhs and Weib (1969). The shallow foundation 
had a length of 2 m and a width of 0.5 m. The water table was at the ground 
surface, the submerged unit weight of the sand was 0.95 t/m 3 , the friction 
angle measured by the triaxial test was 40° to 42°, and the base of the footing 
was 0.5 m below the ground surface. In the first test the foundation was loaded 
centrally and vertically, and failure occurred as a load of 190 t. Hansen sug¬ 
gested the use of </>„. of 40° in his discussion of the Muhs-Weib tests. 

<!> P i = 1-1(40) = 44° 

Eq. 7.4: N q = (20.78)(5.55) = 115.3 

Eq. 7.6: N y = (1.5)(114.3)(0.9657) = 165.6 

Using the first two terms of the Eq. 7.33, failure along the short sides, with 
no inclined loading, inclined soil surface, or cohesion, and noting that d = 1. 
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^ = \ yN y Bs yB + -yD f N q d qB s qB 

Eq. 7.16: s yB = 1 - 0.4 ^ = 0.90 

Eq. 7.17: s qB = 1 + (0.695) ^ = 1.174 

Eq. 7.29: d qB = 1 + 2(0.966)(0.0932)(0.785) = 1.141 

Q d = 1^(0.95)(0.5)(165.6)(0.90)(1) + (0.95)(0.5)(115.3)(1.174)(1.141)J 
= 35.4 + 73.4 = 108.8? 

The second test was performed with a foundation of the same size, in the 
same soil and position as the water table. The foundation was loaded centrally 
to failure with a vertical component V = 108 t and a horizontal component 
in the L direction H L = 39 t. The revised Eq. 7.7 for the loading indicated is 

= 2 yLN y s yL d yI i yL + yD f N q Sq L dq L iq L \ 


because H B = 0; d qB = 1.414; s qB = 1.174; i yB = i qB = 1; Bi yB = 0.5; then 
d qB s qB i qB = (1.141X1.174)(1) = 1.340. 

In the L direction, the following quantities are found: 


Eq. 7.29: 
Eq. 7.11: 

Eq. 7.10: 



Because Bi yB > Li yL , use Eq. 7.33 and the following computations. 
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Eq. 7.23: 

V 

Eq. 7.21: 

S qL 


^qL S ql}qL 


dqB S qBiqB 


1 - 0.49 


1 + 0.695 


0.466 

0.500 

0.739 


= 0.627 


= 2.027 


v 0.500y 

(1.044)(2.027)(0.3695) = 0.782 
(1.141)(1.174)(1) = 1.340 


Because 0.782 < 1.340, use the lower g-term in Eq. 7.33. 


Q d 


1 (0.95)(165.6)(2)(1)(0.627)(0.233) 

+ (0.95X0.5X115.3X1.044X2.027X0.3695) 
23.0 + 42.8 = 65.8 t 


In the hrst test noted above, the computed failure load was 108.8 t versus 
an experimental load of 180 t, yielding a factor of safety of 1.65 based on 
theory alone. In the second test, the computed failure load was 65.8 t versus 
an experimental load of 108 t, yielding a factor of safety of 1.64 based on 
theory alone. 

Of interest is that Hansen employed a value of $ of 47°, instead of the 44° 
suggested by his equations, and obtained excellent agreement between the 
computed and experimental values of the failure load. Muhs and Weib (1969), 
as noted earlier, reported the measured friction angle to be between 40° and 
42°. Had Hansen elected to use 42° instead of 40° for <^> fr , the value of <j> , for 
use in the analytical computations would have been (42)(1.1) or 46.2 degrees, 
fairly close to the value that Hansen found to yield good agreement between 
experiment and analysis. 

Ingra and Baecher (1983) discussed uncertainty in the bearing capacity of 
sands and concluded: “When the friction angle is imprecisely known (e.g., 
having a standard deviation greater than 1°), the effect of uncertainty in 4> 
predominates other sources.” 

Examples from Selig and McKee (1961) In contrast to the footings tested 
by Muhs and Weib (1969), Selig and McKee (1961) tested footings that 
ranged in size from 2 by 2 in. to 3 by 21 in. The soil employed was an Ottawa 
sand that was carefully placed to achieve a uniformity of 112.3 lb/ft 3 in a 
box that was 48 in. square by 36 in. deep. The friction angle of the sand was 
measured in the triaxial apparatus and was found to range from 38° to 41°. 
An average value of the friction angle was taken as 39.5; thus, using the 
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Hansen recommendation, the following value was computed for use in the 
analysis: 


4* =11 ^ = 43.45° 

Some small footings were tested with no embedment, which is not consis¬ 
tent with practice. The experimental results greatly exceeded the results from 
analysis using the Hansen equations. Selig and McKee performed three tests 
where the base of the footings was placed below the ground surface. The 
following equation is solved in the table below, with the following constant 
values: N q = 106.0; N y = 149.2; y = 0.0650 lb/in. 3 ; s y = 0.6; d y = 1.0; s q 
= 1.688; and d q is as shown in the table. 

X = \ + y D f N <, s A, 

For the three experiments that were evaluated, the factors of safety in the 
equation that were computed in Table 7.2 were 1.50, 1.38, and 1.23, yielding 
an average value for a fairly small sample of 1.37. The selection of a value 
of 4> P i is important, and Hansen’s correlation with </>„. may be used unless a 
more precise correlation becomes available. A factor of safety appears in the 
Hansen procedure, as noted from the small sample, and additional analyses 
of experimental data of good quality, particularly with full-sized footings, will 
be useful. 

The settlement at which failure is assumed to occur is of interest when 
performing experiments. With regard to the bearing capacity of the base of 
drilled shafts in sand, several experiments showed that the load kept increas¬ 
ing with increasing settlement, and it was decided that the ultimate experi¬ 
mental load could be taken at a settlement of 10% of the base diameter. A 
possible reason for the shape of the load-settlement curve was that the sand 
kept increasing in density with increasing settlement. 

Muhs and Weib (1969) apparently adopted 25 mm as the settlement at 
which failure was defined. Only a few of the load-settlement curves were 
reported by Selig and McKee (1961), but those that were given showed failure 
to occur at a relatively small settlement, with the load remaining constant or 
decreasing after the maximum value was achieved. In both sets of experiments 
discussed above, the relative density of the sand was quite high, probably 


TABLE 7.2 Computation of Q d to Compare with Experimental Value 


Footing size, 
in. 

Area, 

in. 2 

D f , 

in. 

d . 

Comp. Q d , 
lb 

Exp. Q d , 

lb 

Factor of Safety 
in Equation 

3x3 

9 

2 

1.123 

261 

391 

1.50 

3x3 

9 

4 

1.171 

516 

711 

1.38 

3x3 

9 

6 

1.205 

782 

964 

1.23 
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leading to the maximum load at a relatively small settlement and a settlement¬ 
softening curve for greater settlement. 

The small number of experiments where full-sized footings were loaded 
to failure is likely due to the cost of performing the tests. Furthermore, the 
cost of footings for a structure is likely to be small with respect to the cost 
of the entire structure, so the financial benefits of load tests would likely be 
small in most cases. 


7.5 EQUATIONS FOR COMPUTING CONSOLIDATION SETTLEMENT 
OF SHALLOW FOUNDATIONS ON SATURATED CLAYS 

7.5.1 Introduction 

As noted earlier, the geotechnical design of shallow foundations is presented 
in Chapter 9. Two kinds of settlement are discussed: short-term or immediate 
settlement when the settlement of granular soils (sands) and cohesive soils 
(clays) will occur rapidly, and long-term settlement of clays that can occur 
over a period of years. Some long-term settlement of sands can occur. The 
most evident case is when the sands are subjected to vibration, causing the 
grains to move and increase the density of the soil. Large amounts of settle¬ 
ment can occur if very loose sands are subjected to vibration. Loose sands 
below the water table cannot rapidly increase in density when affected by an 
earthquake and are subject to liquefaction because the water cannot imme¬ 
diately flow out. When a stratum of sand is subjected to sustained loading, 
some settlement of the foundation may occur that is time dependent. 

Long-term settlement of clay deposits below the water table is of concern 
to the geotechnical engineer. Two types of predictions must be made: the total 
settlement that can occur over a period of time and the time rate of settlement. 
The procedures available for making the two predictions are presented in the 
following sections. 

The assumptions made to develop the equations for computing the total 
settlement of a stratum of saturated clay are as follows: 

1. The voids in the clay are completely filled with water. 

2. Both the water and the solids of the soil are incompressible. 

3. Darcy’s law is valid (see Chapter 3); the coefficient of permeability k 
is a constant. 

4. The clay is confined laterally. 

5. The effective and total normal stresses are the same for every point in 
any horizontal plane through the clay stratum during the process of 
consolidation (Terzaghi, 1943). 

Further, Terzaghi noted that the process of consolidation, as water flows 
vertically though soil, is analogous to the heat flow from an infinite slab if 


Copyright © 2006 John Wiley & Sons Retrieved from: www.knovel.com 



214 THEORIES OF BEARING CAPACITY AND SETTLEMENT 


the slab is at a constant temperature and the boundaries are maintained at a 
different temperature. 

If the assumptions noted above are correct, the equations for total settle¬ 
ment and the time rate of settlement of a stratum of clay below the water 
table when subjected to load are perfectly valid. The assumptions made in 
developing the theory cannot all be satisfied; for example, the loaded area 
must be finite, and some water could flow from the stratum laterally rather 
than vertically. However, if care is used in sampling and testing the stratum 
of clay, the consolidation equations can be used, as evidenced by past expe¬ 
rience, to make predictions for the total settlement with reasonable accuracy 
and prediction of the time rate of settlement with less accuracy. 


7.5.2 Prediction of Total Settlement Due to Loading of Clay 
Below the Water Table 

Figure 7.9a shows a stratum of sand below the ground surface above a stratum 
of clay. The water table is shown in the stratum of sand. The assumption is 
made that the clay has always been saturated. The line labeled <Xq shows the 
distribution of pressure from the self-weight of the soil prior to the imposition 
of a surface load, and the quantity Act' shows the increase in cr' due to the 
surface loading. The curve in Figure 7.9b shows the results from a laboratory 
consolidation test of a sample of superior quality from the stratum of clay. 
The dashed line shows the direct results from the laboratory data, with the 
solid line showing the corrected laboratory curve to yield a field curve, as 
described in Chapter 3. The portion of the curves marked “rebound” indicates 
that the soil will swell when load is removed. The simple model for com¬ 
puting total settlement is shown in Figure 7.9c. 

The determination of the change in void ratio Ae for use in computing 
settlement S is shown in Figure 7.9b The average value of Act' shown in 
Figure 7.9a is added to cr f j, and the value of Ae is obtained as shown. The 
total settlement can now be computed with the equation in Figure 7.9c. 

The equation for settlement S in Figure 7.9c may be rewritten in a con¬ 
venient way to compute settlement based on the plot in Figure 7.10. The 
curve in the figure has been corrected to offset remolding and other distur¬ 
bance. The plot of void ratio e versus the log of effective stress cr' is for a 
sample of clay that had been preloaded. The first branch of the curve indicates 
the reloading of the sample after removal of previous stress, the reloading 
curve. The other branch represents the behavior of the sample after exceeding 
the magnitude of the previous loading, the virgin compression curve. The 
equation for the two branches of the curve is as follows: 


5 = 


H 


1 + e n 


log -j + 


(Tq + Act' 


(7.36) 
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Surface Load 


Sand 

Clay 

Sand 



(c) 

Figure 7.9 Data required for computing total settlement of a clay stratum due to an 
imposed surface load, (a) Solid profile showing initial effective stress a, J 0 , surface load, 
and distribution of imposed load with depth A cr'. (b) Field curve of void ratio e versus 
effective normal stress cr'. (c) Model for computing total settlement. 
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log a’ 



Figure 7.10 Plot of void ratio versus log of effective stress for an overconsolidated 
sample. 


where 


c r = compression index (reloading), and 

c c = compression index (virgin loading), dimensionless. 


In some instances, the settlement due to the reloading portion of the curves 
is deemed to be so small as to be negligible. For normally consolidated clay, 
only the virgin compression curve will be used. 

To obtain the quantity A cr', a method for predicting the increase in load as 
a function of depth must be developed. An appropriate method has been 
proposed by Boussinesq (1885) as shown in the following equation (Terzaghi 
and Peck, 1948): 



where 


a,. = stress in the vertical direction, 

Q = concentrated vertical load, 

Z = vertical distance between a point N within an semi-infinite mass that 
is elastic, homogeneous, and isotropic, and 
r = horizontal distance from point N to the line of action of the load. 

The assumption is made that the applied load at the ground surface is perfectly 
flexible, with a unit load of q. The loaded area is divided into small parts, as 
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expressed in the following equation, where load dq is assumed to be acting 
at the centroid of the area dA\ 


dq = qdA 


(7.38) 


Substitution can be made into Equation 7.36 with the following result: 


da v 


Jq_ 

2ttz 2 


1 

1 + ( r/z) 2 


2 

dA 


(7.39) 


The magnitude of a v at point N due the entire loaded area can be found 
by integrating Eq. 7.39. Newmark (1942) produced a chart for easy use in 
the determination of cr v (Figure 7.11). The following procedure is employed 
in using the chart to find cr v at depth z: the distance to the depth is given by 
the length A-B in the chart; employing the scale, a sketch is made of the 
loaded area; the point in the sketch under which the value of cr v is desired is 
placed over the center of the sketch; the number of subdivisions n, including 
fractions, covered by the sketch is counted; and cr v = 0.005 n <x 0 , where <x 0 is 
the vertical stress at the ground surface due to the applied load. The number 
of subdivisions in the chart is 200; therefore, if all of the subdivisions are 
covered, cr v = cr 0 . 

Sketches of the loaded area are made for increasing depths using the A-B 
scale, with the sketches of the loaded area becoming smaller and smaller. 
Thus, a plot of a v as a function of depth z beneath the ground surface can 
readily be made. As shown in Figure 7.12, the value of cr„ as a function of 
depth z can readily be made for a point under a particular footing if a number 
of closely spaced areas exist with differing values of cr 0 . The scale A-B was 
set at 8 ft; thus, the value of <x„ was to be determined at a depth of 8 ft below 
the base of the footings. As shown in Figure 7.12, the scale A-B was used 
to draw sketches of the footings 6 by 6 ft and 8 by 8 ft, separated by a 
distance of 3 ft. The value of cr v was desired below the center of the smaller 
footing. As shown in the figure, the unit load on the smaller footing was 2 
k/ft 2 and the load on the larger footing was 3 k/ft 2 . The number of values of 
n under the smaller footing was 42.4, and n for the larger footing was 10.4. 
The value of cr v can now be computed: 

(cr„)S, = (0.005)[(42.4)(2.0) + (10.4)(3.0)] = 0.424 + 0.156 = 0.58 k/ft 2 

Examining the plotting of the footing sizes in the figure shows that there 
was a lack of precision that had only a minor effect on the result. The New¬ 
mark chart can be used to obtain the increase in effective stress as a function 
of depth beneath a given point under a loaded area due to various horizontal 
spacings of footings, considering the surface loading from each of them. 
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Figure 7.11 Newmark chart for computing vertical stress beneath a loaded area. 
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7.5.3 Prediction of Time Rate of Settlement Due to Loading of Clay 
Below the Water Table 

Terzaghi (1943) reasoned that prediction of the time at which a certain amount 
of settlement occurs could be simplified by assuming that the surface loading 
occurred instantaneously, that excess porewater pressure would develop 
throughout the stratum of clay, and that drainage from the clay would occur 
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instantly where the clay interfaced with a porous stratum of sand. Thus, water 
in the stratum of clay in Figure 7.9 would drain upward from the top of the 
clay and downward from the bottom of the clay. Terzaghi noted that the 
process was precisely analogous to the heat-flow problem, where a semi¬ 
infinite slab of material at an elevated temperature is suddenly subjected to 
zero temperature, or to a lower temperature, at both surfaces of the slab. 

The differential equation for heat flow is shown in Eq. 7.40 (Carslaw and 
Jaeger, (1947) 


du d 2 ll 

dt dz 2 


(7.40) 


where 

u = temperature (excess porewater pressure), and 

k = diffusivity reflecting the ability of the slab of material to transmit heat. 

The problem is to find the properties of the soil, plainly involving the coef¬ 
ficient of permeability k, that will replace the diffusivity k. The term k is 
replaced by c v in consolidation theory, where c„ is called the coefficient of 
consolidation, as defined by the following equation: 


k 

™ v y w 


(7.41) 


where the rate of flow of water from the clay is directly proportional to the 
coefficient of permeability k and inversely proportional to y w and m,„ where 
m v = aj\ + e and a v = —de/dcr'. Thus, if the coefficient of compressibility, 
a v , is larger, meaning that the coefficient of volume compressibility, m v , is 
larger, the rate of flow will be less. Thus, the differential equation for the rate 


at which water will flow from a stratum of clay that has been subjected to an 
increase in effective stress, cr', due to the placement of a surface loading 
above the clay, is shown in Eq. 7.42: 


dll d 2 ll 

dt ~ ° v dt 2 


(7.42) 


The solution of Eq. 7.41 requires the implementation of the following 
boundary conditions: the excess porewater pressure, u, is equal to zero at the 
drainage faces; the hydraulic gradient i is equal to zero at the mid-height of 
the stratum of clay and at the boundary where flow is prevented, where i = 
duldz', and after a very long time (infinite according to the differential equa¬ 
tion), u = 0 at all depths. 
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The solution of Eq. 7.42 for the conditions given is as follows: 

U z (%) = f{T u , £) (7.43) 

where 

t (7-44) 

The dimensionless number T v is called the time factor. The solution of Eq. 
7.43 in terms of dimensionless coefficients is shown in Figure 7.13. The figure 
illustrates the process of consolidation graphically by showing the distribution 
of excess porewater pressure, u, with the passage of time. The figure readily 
shows that a very long time will be required for u to equal zero throughout 
the stratum. The hydraulic gradient causing the discharge of porewater be¬ 
comes smaller and smaller over time. 

The data in Figure 7.13 can be integrated to produce Figure 7.14. The time 
for a given percentage of the consolidation to take place in the laboratory for 
a representative sample is determined, usually at 50% of the laboratory con¬ 
solidation. Using the time factor T v of 0.197 (see Figure 7.14), the value of 
the coefficient of consolidation c v can be found. And Eq. 7.43 can be used 
to plot a curve showing the amount of settlement of the clay stratum as a 
function of time. 



Figure 7.13 Percentage of consolidation U z (%) as a function of relative depth Z/H 
and time factor T v (note: thickness is 2H for the drainage top and bottom of the 
stratum) (from Peck et al., 1974). 
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Time factor, T Y (log scale) 


Figure 7.14 Theoretical relationship between degree of consolidation U and time 
factor T v (from Peck et al., 1974). 


PROBLEMS 

P7.1. a. Compute the total load, Q q , in kips, using the Hansen equations, for 
a footing 4 ft wide by 8 ft long founded on a stratum of sand. The 
vertical load is uniformly distributed; there is no horizontal load. 
The value of D f is 3.5 ft, the water table is at the base of the footing, 
and the friction angle from the triaxial test is 30°. For the sand above 
the base of the footing, the water content is 20%, and the unit weight 
y is 120 pcf. Compute a value of submerged unit weight y' if 
needed, assuming specific gravity of the particles of sand to be 2.67. 
b. The problem remains the same as in part a except that the footing 
rests on overconsolidated clay at a depth of 3.5 ft and the footing 
rests on the clay. Assume the same unit weights as before and de¬ 
termine the undrained shear strength of the clay to yield the value 
of Q q found in part a. 

P7.2. a. Under what condition might you wish to construct a footing with a 
base at some angle with the horizontal? Show a sketch, 
b. What technique could you use to fine the most favorable angle for 
the base? 

P7.3. Under what conditions would you wish to use the Hansen equations 
for a sloping ground? Show a sketch. 

P7.4. The time to reach 55% consolidation for a laboratory sample of clay 
that was 0.6 in. thick and tested under double drainage was 24.5 sec¬ 
onds. How long would it take for a stratum of the same soil in the held 
that was 14 ft thick and drained on only one side to reach the same 
degree of consolidation? 
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